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Abstract. We consider the inhomogeneous Dirichlet problem on product do- 
mains. The main result is the asymptotic expansion of the solution in terms 
of increasing smoothness up to the boundary. In particular, we show the exact 
nature of the singularities of the solution at singularities of the boundary by 
constructing singular functions which make up an asymptotic expansion of the 
solution. 



1. Introduction 

In [2], we examined the <9-Neumann problem for (0, l)-forms on product domains 
in C™ of the form D = D± x • • • x D n where C C. In our analysis we related 
the solution of the 9-Neumann problem to the inhomogeneous Dirichlet problem, 
and as a corollary obtained specific information on the nature of the singularities. 
We were able to construct singular functions which we used to write an asymptotic 
expansion of the solution, each successive term in the expansion exhibiting a higher 
degree of differentiability up to the boundary. This paper is written to generalize 
this corollary result of [2]- Although we focus on the inhomogeneous Dirichlet 
problem, the methods used can be applied to other boundary value problems as 
well. 

We let f2 C R™ be a product of q smooth bounded domains, fi = Oi X • • • X f2 g , 
where £li C R +J< , for 1 < % and = n — q. Locally, in a neighborhood, U 

of xq, a point in the distinguished boundary of ft, dfti x • • • x dfl q , pi is a defining 
function for fij, which we shall assume to be of the form 

Pi =<t> i (t{,...,t) i )-x i 

for 1 < i < q, and fa S C°°(R j >). Thus, Q n U is the set of all x e fl such that 
Pi(x) < for all 1 < i < q. 

We consider the inhomogeneous Dirichlet problem on the product domain, Cl 

(1.1) Au = f inQ 

u = on dil 

and the singular behavior of it in a neighborhood, U , of a point, xo, of the distin- 
guished boundary, at which dfl is not smooth. In our analysis / will be in the class 
of C°°(0). Let H a (il) denote the Sobolev a space. Existence and uniqueness of a 
solution follows from Jerison-Kenig 0]: 

Theorem 1.1 (Jerison-Kenig). Let Q be a bounded Lipschitz domain in W 1 . Sup- 
pose that j < ol < f ■ Then the inhomogeneous Dirichlet problem U.l)) has a unique 
solution u <E H a {Q) and 

\W\\H"(Q.) ^ Il/llff°-2(f2) 
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for every f 6 H a 2 (Sl). 

In the same paper 0] (see also [5]) Jerison and Kenig also proved that in our 
situation, u E (Q). In contrast to the case of smooth boundary, the non-smooth 
case exhibits singularities. Thus the classical L 2 treatment (see Lions-Magenes [B]) 
in which gains in derivatives are obtained has to be modified. In this paper it is 
our purpose to write an explicit solution as a sum of terms with increasing degrees 
of differentiability up to the boundary, and thus give an analysis of the behavior of 
the solution near the singular parts of the boundary. 

In [5] , the case of ttj € R 2 for each j was considered and a conformal mapping 
was used to reduce the problem to one on the product of half-planes. Instead of a 
conformal mapping, in this paper we use a change of coordinates, so that, locally, 
the domain looks like a product of half-spaces. 

We organize our paper as follows. In Section 2, we setup the problem, and trans- 
form it to one on a product of half-planes. We begin the process of localizing the 
problem, so that we only concern ourselves with function behavior in a neighbor- 
hood of a singular point of the boundary. In Section 3 we write down an infinite 
sum, which represents a solution modulo terms smooth up to the boundary for the 
problem setup in Section 2. In Section 4 we construct explicit singular functions 
which comprise the terms in the asymptotic expansion of our solution. 

2. A PROBLEM ON A PRODUCT OF HALF-PLANES 

Recall from above that ft = fii x • • ■ x il q , and in a neighborhood of a point 
on the distinguished boundary has as a defining function 

4>i \t\ > • • • i tj t ) 2-i ■ 

Wc thus use a transformation of coordinates 
(2.1) y =Xj -<f>j 1 < j < q 

Vq+ra =t l k 1 < m < n - q, j\ H h ji-i + k = m, k < ji 

where xq corresponds to y = 0, and the related matrix A = [dij], where 

f 1 i = j 

aij = < -4>j,k i>q,j = i 

and 4>j t k — where k and / are determined by the unique representation of i as 

i = ji + . . . + + k, and k < j t . 

The transformation ()2.1ll leads to a Dirichlet problem on the domain H™ in W 1 
given by ij = i 1 x ■■■ x i« x where W = {(y u . . . , y n ) : y { > 0} for 

i = 1, . . . , q in which the differential operator, A, is replaced with 

n d 2 ™ " dak d 

£, 9 d Vl dy 3 + § 5 Qhl % d yj 
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where g lJ = Y^k=i a ki<ikj is the metric tensor given by the entries of the matrix 
A 1 A. Because of the invertibility of A, and from the smoothness up to the boundary 
of the we see that € C°°(H?) as are the terms dki^S^-- Calculating g % i gives 



A'=E(1 + |V^:| 2 )^ 



1=1 



i=l j,k=l 



2 EE ( wj ()!hrj . hkdy . + EE*i d!h lh/ 



EC . \ ^ \ - oa ki O 
7^2+ 1^ 1^ ak r 



Therefore, using the change of coordinates (|2.ip we examine the problem 
tpA'u = ipf in H™ 

u = on 2/j = 0, i = 1, . . . , q, 

where ip € C^°(il) is cutoff function such that < <p < 1 and supp ip C U, and 
93 = 1 near ieo- 

Commuting ip with the operator A' gives 

(2.2) AVu = <pf + [A', 

We begin our discussion of the singularities of the solution by noting the location 
of the singularities must be along intersections of boundaries of more than one d£li. 

Lemma 2.1. Let f S C°°(Q), u S L 2 (fl) be the unique solution to the Dirichlet 
problem, II 1.1]) . on fl. Let V C such that V D dilj 7^ /or at most one j . Then 

Proof. The lemma is a consequence of regularity at the boundary of the Dirichlet 
problem [3]. □ 

3. An asymptotic construction of a solution 

We set v = ipu and h to be the right hand side of (|2.2p . and note that h = f 
in a neighborhood of xo- Without loss of generality we can assume h is just ipf 
since the error term [A', ip] u leads to terms smooth up to the boundary in a small 
neighborhood of our chosen point xq. This can be seen by writing the solution 
in terms of Green's function. Let 4>{y) be a function which is equivalently near 
xq, and let G(x,y) be the Green's function for the Dirichlet problem on S7, then 
G(x, y)<j)(y)dV \y) can be extended smoothly across the boundary near Xq- 
We shall make use of odd reflections along y.- L = for 1 < i < q. We take odd 
reflections of (|2.2p with respect to the variables yi for 1 < i < q and we denote the 
extension by a superscript o. We also denote by a superscript Oj an even extension 
in variables yi for 1 < i < q, i ^ j and an odd extension in yj, and a superscript e 
is used to denote even extensions in all variables yi for 1 < i < q. We write 1|2.2|) . 
after reflections, in the form 
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Due to our discussion above on the error term in (12. 2| leading to terms which are 
C°° smooth up to the boundary of SI in a neighborhood of xq, we look to solve 
(j3~Tj) with h° replaced by (ipf)°. 

We make the note here that the difficulty in following the process in the case 
of smooth domains to attempt to construct a parametrix lies in the existence of 
singularities in the symbol along yi = for 1 < i < q due to the reflections. This 
difficulty can be resolved by referring to Lemma 12. 1( and we could then proceed 
to show the error terms of the parametrix construction yield an error term in 
the solution which is also smooth up to the boundary as we do in Theorem 13.21 
below. Since this examination of the error terms is independent of the method 
of parametrix construction, we choose here to follow our analysis in [2] from the 
beginning, which nonetheless has the flavor of a parametrix construction. 

Here rji will be the Fourier variable corresponding to yi for i < q, and £j the 
Fourier variable corresponding to yi for q + 1 < i < n. We denote the symbols of 
bi and c, by Bi and Ci, respectively, for 1 < i < q. Let (q, i) = q + ji + ■ ■ ■ + ji-i 
and from the symbol of the operator in (|3.ip . 



i i i 

P(v,t,v) = - E( 1+ I V< ^I 2 H 2+ E B Tivi+ii ■ ■ ■ ,y n ,0vi-^2 Ci(y q+1 , . . . ,y n ,0, 

i—1 i—1 i—1 

we define the operator K by 

i , 

Ku = -Y^ I (B^rji - d) ue iy <t>rid£dr), 

i=l J 

and with 



where ~ refers to a partial Fourier transform in the variables y\, . . . ,y q , we induc- 
tively define 



f33) t - , - 1 [ y(r,\ Kv 3^Vl+U--^yn) «„!,■■■,„„)■„, >Q 



We write the solution v° to (|3.ip in the form 

(3.4) v° = v H h v N + v Rn 

for any integer N > where the remainder term, vr n , satisfies 

9yf 1 dyi 




Kv N + F.T.- 1 ^(1 - x) 

We write F.T.^ 1 to refer to an inverse Fourier transform in the rj variables only. 
The forms of the Vj are written in the following lemma. 

Lemma 3.1. Each Vj is of the form 



fSfiA ,,.-/ Yi+^ri) i2^S±llllll^l5±l> (<of) eiv^dfdn 
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where Pj{y, £, rj) is C°° smooth as a function of y q +\, . . . , y n and is a polynomial in 
£ and rj, of order 5j in the n variables. 

Proof. The lemma would be clear if it were not for the odd reflections in the sym- 
bols, B k . We thus examine terms of the form H (y k )rj k u° , where H(t) is defined to 
be 1 for t > and -1 for t < 0. We have 

F.T.- 1 {n k u°)^ l — + tS(y k )u°, 

oyk 

where e k denotes an even reflection in the y k variable, and odd reflections in 
2/1, ... , Vk-uVk+i, •■■iVq- And thus 

H(y k )F.T.^ ( Vk u°) = i + i%*K. 

Using 



I I = —irj k u eh 



( du 
\dyk. 

in (|3.3p . and noting that, since Vj is odd in all r\ variables by our construction, 
v j(yk = 0) = for k = 1, . . . , q, we can write 



"X^ 1 + l v ^| 2 )^^+i = \ B k(yg+U ■■■,y n ,0 (iVkVj" +iF.T. y - k Vj(y k = 0)) 
fc=i fe=i 

+ C k (y q+1 , . . . ,2M,£)£j) 

(3-7) = ^ (£ fe (?/ g+ i, . . . ,y n ,O ir )kVj k + C k {y q+1 , . . .,y n ,Qv. 

fe=i 

where F.T. VI refers to a Fourier transform in all y\, . . . , y q variables but y k - Now 

dyi 



+ \W4>k\ 2 )ril{-im^ 1 ) = -imxWKvj-* +i m F.T. y - h ^{yi = 0), 



k=l 

and so 

1 J J " ELi(i + NM 2 H + ELi(i + |v<fe| 2 )^ 3 ~ L > 

while for vq we have 



(3.9) 



ELi(i + |v^| 2 )^ ELi(i + lv^l 2 )^ 



wo, from p.2[) . has the form of (|3.6p and taking an inverse transform with respect 
to rji of the last term on the right side of (|3.9[) shows that it also has a Fourier 
transform in y\ , . . . , y q of the form in (|3.6p . We use an induction argument and 



assume vj is of the form in (13.61) 



The proof will be complete, by (|3.7p . if we show v^ k is also of the form (|3.6p . 
An inverse transform with respect to n k of the last term on the right side of (|3.8p 
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shows that it also of the form in (13. 6|) . For the first term on the right of (|3.8| we 
use 

------ e, q r 

Kvj-x = - RT -v J ( B Tm - C,) vf^y^drjd^ 



3-1 



1 , 

= -J2 F.T. n J {B m - d) vfl^-^dr,^ 



- (Birji - Ci)v 3 -i, 

where en means an even extension in both yi and yi, and we use the fact that 
an even extension in yi of a function H(yi)F.T.^ 1 (rjiu°) restricted to yi > has 
a transform rjiu° if u{yi = 0) = 0. The same arguments above showing v^ k is of 
the form (|3.6[) shows we may apply the induction hypothesis to Vj l lx as well. The 
induction hypothesis then shows that v^ also has the form in (|3.6[) . and with Vj l 
inserted into ()3.7|) we have Vj+\ of the form in ()3.6|) . □ 

We show that (|3.4[) gives an asymptotic expansion of our solution v. 

Theorem 3.2. Each successive term in \3.4\ ) is of increasing differentiability up to 
the boundary, and the remainder terms Vr n are also of increasing differentiability 
up to the boundary. The representation \3.4\ ) is therefore an asymptotic expansion 
of the solution to the Dirichlet problem. 

Proof. From our definition of the Vj , an induction argument shows each Vj is infin- 
itely differentiable with respect to the variables y q +i, Furthermore, it fol- 
lows by our construction of the Vj, and by (|3.6p . that if Vj € H k (M. n ), the Sobolev-fc 
space, then Vj+i £ H k+1 (R") . The proof of Lemma |3~T1 also shows that restricting 
Vj to yi > for 1 < i < q and then reflecting about yi = in an even manner to 
form gives v^ 6 H k (M. n ). To see the remainder terms, vr n are of increasing 
smoothness up to the boundary with increasing N , we take Fourier transforms of 
(13. 5[) with respect to the variables y\ , . . . , y q to obtain 



N-l 




(3.10) F.T.- 1 (P(y q+1 ,... 1 y n ,^ V )v RN ) = Kv N + (1 - X ) 

From the proof of Lemma 13. 1[ we have 
-— — - q 

Kv N = ^ [B k (y q+l , . . . ,y n ,QirikVN + C k (y q+ x, ■ ■ ■ , y n , O^n) ■ 



fc=i 



If we assume vn £ 7J fc ( Ar ^(K"), and hence v^ £ H k ( N ^ (W 1 ) from above, for some 
k(N) strictly increasing in N, with the property that ^wjv £ iJ fe ( JV ) (R") for any 
multi-index a — (at, . . . ,a n with a\ = • ■ ■ = a q = 0, then Kvn £ H k ( N ' > ~ 1 (W 1 ) by 
definition of the operator K, and Pk(N)-i(v)KvN £ L 2 for Pk^-iiv) a polynomial 
in the n variables of degree k(N) — 1. Also, since any polynomial of any degree in 
the r\ variables multiplied by the last term on the right hand side of (|3.10[) is in L? 
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due to the compact support in 77 of the term 1 — x, we see, by multiplying (|3. 10|) 

by Pk(N)-i(v), that 

d d 



is the solution to a problem with mixed Dirichlet and Neumann conditions, and 
by 0], its restriction to H™ is in H a {W^) for a < 3/2. Therefore, Uh n |h» € 
# fc W+<*(H£) for a < 1/2, and so is of increasing smoothness up to the boundary 
for increasing N. □ 

4. An explicit calculation of the singularities 

In order to determine the singularities in each term in the expansion (|3.4p . we 
take inverse Fourier transforms of the expression (|3.6[) with respect to the r\ vari- 
ables. 

Without loss of generality we suppose the cutoff function x is °f the form 

Xm ' ' ' Xr/g , 

where Xm is a cutoff function in the variable r\i only, equal to in a neighborhood 
of t]{ = 0. We can do this because 

X 3+1 (v) - Xm ■■■Xn q 

can be written as a sum of terms, each of which has support contained in large rjk 
for at most one k. Hence, 

where a* = 1 + |V<^i| 2 , is a sum of terms which have infinite decay in all but one 
rjj, and thus is the transform of a function which is C°° in all variables but one 
Uj. As such a term is odd in jjj and as the denominator Q3i=i a i r iiY'' +l is the 
symbol of an elliptic operator, such a term is the solution to a Dirichlet problem 
on a half-plane, and is therefore C°° smooth up to the boundary of the half-plane. 
In summary the difference between (|3.6[) and 



(4 1) y ...y MMllllll^lM ^n 

(4.1) Xm X Vq W) 

is the transform of a term which, when restricted to H™, is in C°°(H^). 

We then take (|4.ip and integrate by parts with respect to the variables yi for 
1 < i < q in the Fourier integral of ipf (r],£), starting with y\ and proceeding to y q , 
keeping as remainder terms those which have decay in one rj variable to the order 
—2(N + q) — 1. Such terms are of the form 

(A \ v v Pj(Vq+l,---,V*,Z,V) 1 1 

K^-^i Am ' ' ' ai), " 



/ d 2(N+ g ) + l d a {<pf)°. \ 

F - T -k+i,...,n y d ^ 2{N+q)+1 dya (o,...,o,y fe+1 ,...,t/„) j 

in which a is a q index for which ak+i = ■ • • = ot n = 0. We shall show below 
in Theorem 14.21 that such remainder terms are sufficiently continuous up to the 
boundary of H" . 
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Setting aside the remainder terms (|4.2j) . we analyze those other terms which 
result from the expansion of (iff) (r], £), and we are led to study terms of the form 

where 4>k,j{y, v) takes the form of odd reflections along j/j = for each i = 1 . . . <j> 
of functions of y smooth up to the boundary, and is a polynomial in £ and rj of 
order 5j in the £ and rj variables, and a polynomial of degree j in the r] variables.. 
We use the notation r\ k = rfc 1 ■ ■ ■r\ q Zq for k = (fci, . . . , k q ). 

Up to multiplication by a constant the following relation holds for < I < | 

Jim,-- ,y q )-v 



e nyi,-,y q )-v dv 



« 9 (ELi ntiV^ 



Let 



(4.4) = j-, < I < q 

^i— 1 ttj 

uciine $^ 



For Z > | and g even, we define $? to be the unique solution of the form 



Pi (y) log \J2 +P2(y), 

where ^2(2/1 = 0) = 0, pi and P2 are polynomials of degree 2Z — q in the variables 
yi for 1 < i < q, and are C°° smooth with respect to variables y.; for q + 1 < i < n, 
to the equation 

(4.5) 1^ 
For the case q = 2, we take 

<f 1 = -| 1 °sf- + - 

2 \ai a 2 

For q odd we define <I>* for I > | as in (|4.4|1 . the Fourier transform being calcu- 
lated in the sense of distributions. 

In the sense of distributions, we take transforms of (|4. 5[) and write, for q > 2 



d 



which implies 



= -2i(/-l)oi 

(ELr^ 2 ) 

X^*? = -2(i-l)aix w 



for q > 2. For the case 5 = 2, we use pQ to write 



XviXmX^i = XviXm H 2 + s > 

oi??f + a 2 ri2 

where x is a cutoff function such that x = 1 m a neighborhood of the origin as in 
Section [21 and here and below we use s to denote terms which after taking inverse 
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transforms give C°° functions in some neighborhood of the origin in K™. We have 
therefore established the property for / > 

(ELi a ^t ) 

up to multiplication by a constant. 
For k = (k\, . . . ,k q ) we further define 



to 9 rvi rt 



■■•/ ••■/ •••/ *?(^. ..,£)<•••<•••<#•••<. 

Jo Jo Jo JO 
We have the following key property of the which allow us to match these 
functions with the singularities of our infinite sum solution to the Dirichlct problem 
in JO). 



Proposition 4.1. Let <J>^, be defined as above. Let \ € Cq°(W 1 be a smooth cutoff 
function such that x = 1 in a neighborhood of 0. Up to multiplication of a smooth 
function of yi for q + 1 < i < n, has the property 

( 4 - 6 ) Xm ■ ■ ■ Xr, q X®1 k = • ■ • Xu, 4 ' 



Proof. We have 

(-l) |fc| /" <9 |fc| 



-1)1*1 f <9 |fc| -/ 

/ \-7T*"u' ' " *"<'// ' *• 



where the term s in the last line comes from taking derivatives of the cutoff y. By 
definition 

g\k\ 

Therefore, 

f_l)|fc| . — , 

( 4 - 7 ) x*& = ^x*? + ; 

and from the definitions and discussion above, 

= ^ + s 

(4-8) = <Ky ? +i, • • • i 2/«)t^ + s 

(ELi o<»7f ) 

for some smooth function <p. 

Inserting (|4.8|) into (14. 7[) finishes the proof. □ 

For a multi- index p = (pi , . . . , pi) , in which < g, we define $?. for 2 < i < q — 1 
in the same fashion as we did <& q lk but with respect to the i variables y Pl , . . . , y Pi , 
in particular A: is a multi- index of length i. Thus 

( 4 - 9 ) Xr, P xK = X Vp -T , , rj + s, 
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where Xn P = Xn P1 ' ' ' Xv Pi ■ We also use the notation $f fc = when p= (l,...,q). 
From (14. 6|) , (|4.9p , and the construction of the Qf k we see that 

■<* (x^,(xi£)-xi£) 

for large enough m, mi, . . . , m n . Now %jj p — 1 has support for large rjk for at most 
one k, and thus rj" 11 ■ ■ ■ rff"* yXn v (x^ffc) — X^uJ is the transform of a solution to a 
Dirichlet problem on a half-plane with data smooth up to the boundary. Then, by 
the ellipticity of the operator whose symbol is (Ei=i o-ifif) > we obtain derivatives 
of the inverse transform of 

x*(x*&-(x*£) 

is a C°° function on a half-plane. By inverting derivatives, by integration, and 
invoking Lemma 12.11 to show the constants of integrations are smooth up to the 

boundary, we see the terms described by Xv P (.X^fk) an< ^ (x^/k) differ by functions 
smooth up to the boundary, and X&fk wm thus be seen to describe the singularities 
of the solution. 

With a slight abuse of notation we also use the notation §>f k even after a change 
of variables back to the domain il. We are now ready to prove the 

Theorem 4.2. Let f 6 C°°(£l) and let u € L 2 (0) be the unique solution to the 
inhomogeneous Dirichlet problem on fl. Then near the distinguished boundary, 
<9f2i x • • • x d£l q , u is of the form 

(4-10) u= ]T c klp ^ k 

|fc|>0,(>l, 

2<t( P )<q 

where Cki p S C°°(Q), and where are defined as above. 

Proof. (|4.10p is obtained by following the procedure outlined above, matching 
terms, (|4.3[) . in Fourier space to the appropriate function We thus need to 

study the effect <j)k,j has as a polynomial in the £ and rj variables in (|4.3[) on the 
functions <I>^ . We also have to treat the remainder terms (|4.2|) . 



We first show the remainder terms given by (|4.2|) are described by the functions 
<&f fc . We consider the case in which (|4.2j) is given by 



l 



111 / d 2(N+q '> + 

^ ' ' ' Xl N 1 ---%-i ?? f v ^+ 1 ELi«^ F ' T ' m \dyT +q)+M) ° ) 

restricted to yi = for 1 < i < </ and i ^ m, where F.T. m denotes the partial 
Fourier transform with respect to y m . The other such remainder terms are handled 
in a similar manner. The fraction 

1 



is expanded in a geometric series in 

a m r l m 
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up to (N + q)/2 terms (we assume (N + q) is even), the first terms leading to the 
functions <I>f fc in which m ^ p, while the last term is, up to multiplication by a C°° 
smooth function of the variables yi for q + 1 < i < n, given by 
(4.11) 



F.T^Ip^/))' 



Xrn'-'Xv 9 ^N+q+1 , v (iV+g)/2+l 



E?=i a in 



2 ' 



The term in parentheses is the transform of a function odd in the variable y m and 
in C N '(R 9-1 x {y m > 0}), and so all of (|4 . may be viewed as the solution to a 
Dirichlct problem depending C°° smoothly on the parameters yi for q + 1 < i < n 
on a half-space, and by regularity of the Dirichlet problem, the term in (|4.1ip , after 
an inverse transform, is in C (K n ). 

We finish the proof of Theorem l4.2l bv showing that the polynomials of the r\ and 
£ variables in the numerator of (|4.3p still preserve the form of the functions $?, . 
Using the fact that the data function (ipf)° has infinite decay with respect to the £ 
variables, and that multiplication by rji corresponds to differentiating with respect 
to yi the following relations may be used to complete the proof of the theorem. Up 
to multiplication by smooth functions of yi for q + 1 < i < n we have 

d q q 

~Qy~®lk = ®l(ki,...,ki-h-,k q )> 

d 

with the obvious analogies for the Qf k . 

Lastly, we restrict our expansion (|3 .4[) to the product of upper half-spaces, H 9 , 
and obtain an expression of the solution modulo terms in C°°(M. n ) in terms of the 
functions <frf k . After changing coordinates back to f2, we obtain the expansion in 
(00) . ' □ 

Our results here are comparable to those in [2], in which fli C M 2 , and we note 
that an increase in the dimensions of the Qi do not affect the form of singularities 
occurring. See also [T] for the specific case of SI = Sli x fi 2 C R 4 , where fti C M 2 
for i — 1,2, in which logarithmic and arctangent singularities are found along the 
corner. 

Lastly, it is important to mention that there are cases in which the solution to 
the Dirichlet problem does exhibit singularities. The sum (|4.10[) is not trivial; the 
coefficients Chip are not always 0. The example / = 1 on Q reveals this to be the 
case. 
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